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ABSTRACT 



Aims. We interest in the case of a main-sequence star deeply penetrating within the tidal radius of 
a massive black hole. We focus on the compression phase leading to a so-called pancake configu- 
ration of the star at the instant of maximal compression. The aim is to study the tidal compression 
process paying particular attention to the development of shock waves; to deduce reliable esti- 
mates of the thermodynamical quantities involved in the pancake star; and to solve a controversy 
about whether or not thermonuclear reactions can be triggered in the core of a tidally compressed 
star. 

Methods. We have set up a one-dimensional hydrodynamical model well-adapted to the geom- 
etry of the problem. Based on the high-resolution shock-capturing Godunov-type approach, it 
allows to study the compression phase undergone by the star in the direction orthogonal to its 
orbital plane. 

Results. We show the existence of two regimes depending on whether shock waves develop 
before or after the instant of maximal core compression. In both cases we confirm high compres- 
sion and heating factors in the stellar core able to trigger a thermonuclear explosion. Moreover, 
we show that the shock waves carry outwards a brief but very high peak of temperature from the 
centre to the surface of the star We tentatively conclude that the phenomenon could give rise to 
hard electromagnetic radiation, to be compared to some X-ray flares already observed in some 
galactic nuclei harbouring massive black holes. Finally, we estimate that the rate of pancake stars 
should be about 10"' per galaxy per year. If generated in hard X- or y-ray band, several events of 
this kind per year should be detectable within the full observable universe. 

Key words, black hole physics ~ stars: evolution ~ galaxies: nuclei - hydrodynamics - shock 
waves - method: numerical 

1. Introduction 

A massive black hole (BH), of mass can tidally disrupt (small) main- 

sequence stars penetrating within the so-called tidal radius, defined in order of magnitude by 
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with M» and R, respectively the mass and radius of the star With an estimated frequency of one 
stellar disruption every 10"^ years (Magorrian & Tremaine 1 19991 1. or every 10^ years in the pres- 
ence of a self-gravitating accretion disc around the massive BH (Karas & Subr |2007] l. such an 
astrophysical event is expected to provide gas for feeding a moderately active galactic nucleus, and 
for bringing back to active hfe an otherwise quiescent galactic nucleus such as that of the Galaxy. 
Some observations of UV flares (Renzini et al. 11995( 1 and X-ray flares (see Komossa 120021 for a 
review) from the core of non-active galaxies have already been tentatively interpreted in this sense. 

Several analytical and hydrodynamical models have been set up to study particular aspects of 
the tidal disruption process of a main-sequence star by a massive BH. When the star slightly pen- 
etrates within the tidal radius, the dynamics of disruption has been followed e.g. by Nolthenius 
& Katz ( [19821 [T983T l; Regev & Portnoy ( [19871 ); Evans & Kochanek ( [T989l l; Laguna et al. ( [19931) ; 
Khokhlov et al. ( [T993l l; Fulbright ( [T995l l; Diener et al. ([T9971 ); Ayal et al. ( [20001 ); Ivanov & Novikov 
( 120011 ); Ivanov et al. (120031 ); Ivanov & Chernyakova ( 120061) . The distribution and subsequent evo- 
lution of the stellar debris have been followed e.g. by Evans & Kochanek (119891 ); Laguna et al. 
( [19931 ); Kochanek ([19941 ); Lee & Kang ( [19961 ); Fulbright ( [19951 ); Kim et al. ( [19991 ); Ayal et al. 
( 120001 ). A quite complete and coherent picture of the tidal disruption process can be drawn from 
these diflferent studies in this range of moderate star-BH encounters. 

In this paper, we interest in the case when the star deeply penetrates within the tidal radius. 
In such strong encounters, it was noticed on the basis of geometrical and qualitative arguments 
(Carter & Luminet[T982| that before being fully disrupted, the star should pass through a short-lived 
strong compression phase, at the end of which it should adopt a pancake-shape configuration highly 
flattened in its orbital plane. As a consequence, it was suggested that thermonuclear reactions could 
be triggered within the stellar core leading to an explosive disruption. The astrophysical importance 
of such an event would reside in the fact that the additional thermonuclear energy release should 
produce a luminous flare from the core of the galaxy, hence providing an observational signature of 
a pancake star disruption; that it could be sufficient to give the stellar gas a velocity greater than the 
escape velocity from the BH, hence challenging the subsequent accretion scenario (see e.g. Rees 
119881 ); finally that some heavy isotopes could be formed through nucleosynthesis processes which 
would be ultimately ejected in the interstellar medium. 

We shortly review in the following section the main studies on strong star-BH encounters before 
explaining the motivation of the present paper. 

2. The affine star model versus hydrodynamical models 

The first investigations on pancake stars began more than twenty years ago within the framework of 
the semi-analytical affine star model (AM), both in Newtonian dynamics (Carter & Luminet [1983t 
Luminet & Carter [19861 ) and in general relativity with Schwarzschild BHs (Luminet & Marck 
119851 ). The mechanical properties of the AM are described in detail in Carter & Luminet ( 119851 ). 

The main result of the AM was to explicitly show that along a plunging orbit within the tidal 
radius, the tidal field of the BH increases sufficiently rapidly to transitorily surpass the pressure 
and self-gravity fields of the star (Carter & Luminet [1983]). The stellar matter then lands up in 
free fall in the external gravitational field, and compresses in the direction orthogonal to the orbital 
plane in direct response to tidal effects. Finally near the periastron of the orbit, the pressure field 
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suddenly increases within the stellar core to counteract the compressive contribution. The stellar 
matter bounces within the core, and quickly expands in an almost symmetric way to the free fall 
while the star is still moving within the tidal radius. In the relativistic case, the star actually under- 
goes several successive compression phases when its orbit crosses itself inside the tidal radius. 

At the instant of bounce, the density and temperature of the stellar core, initially equal to the 
equilibrium values p, and T,, were found to follow the respective power laws 



considering a polytropic gas of adiabatic index y, and defining the penetration factor within the 
tidal radius 



triggered during the tidal compression was taken into account in great detail in Pichon (|1985l l and 
Luminet & Pichon (fT989] l. 

To go beyond the AM in dealing with the behaviour of the tidally compressed star in a more re- 
alistic way, three-dimensional hydrodynamical models based on smoothed particle hydrodynamics 
(SPH) were developed by Bicknell & Gingold (fT983] l, Laguna et al. (fT993] l, and Fulbright ( fT995] l. 
These models interested in polytropic stars y = 5/3 of solar mass and radius, but did not take into 
account the possible thermonuclear energy generation. The tidal gravitational field was handled in 
Newtonian dynamics, except in Laguna et al. ( 1 19931 1 where it was described in general relativity 
for Schwarzschild BHs. 

The SPH simulations qualitatively confirmed the tidal compression process proposed by the 
AM as a free fall phase followed by a bounce-expansion phase. They nevertheless reported that 
the stellar core could be less compressed than initially predicted, which could possibly impact 
the explosive disruption scenario. In particular, for star-BH encounters with yS < 10, Bicknell & 
Gingold ( I1983I I found that the maximal density of the compressed stellar core scaled as yS'^. Laguna 
et al. (I1993I I confirmed this result, whereas Fulbright (1 19951 1 found a ff"^ dependence. For /? > 10, 
Bicknell & Gingold ( I1983I I found that the maximal core density decreased as p increased. 

The authors of the SPH studies mentioned that the ellipsoidal assumption on which is based 
the AM could partly be responsible for the discrepancy. It is true that the SPH simulations clearly 
showed that the global configuration of the star within the tidal radius was not ellipsoidal. At 
the instant of maximal core compression, the star had more a crescent-like shape within the or- 
bital plane, while it was not uniformly flattened in the orthogonal direction. However, it must be 
pointed out that this configuration essentially concerned the stellar envelope and hardly the stellar 
core. Moreover, since the velocity of the star near the periastron is high enough, of the order of 
(GMbh/^p)'^^ and comparable to the speed of light when the periastron is close to the BH horizon, 
it must also be pointed out that the non-uniform flattening is fully negligible. It is therefore quite 
unlikely that the ellipsoidal assumption of the stellar core postulated in the AM could be rightly 
questioned. 

Noticing an increase of the stellar gas entropy in their simulations, Bicknell & Gingold ( I1983I I 
and Fulbright ( 1 19951 1 also mentioned that shock waves developed during the compression. On the 
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contrary, the AM assumed that the compression process occurs without shock wave development at 
least until the instant of bounce, so that the kinetic energy of compressive motion fully converts into 
internal energy of stellar gas. One could then imagine that the shock wave development prematurely 
halts the free fall phase, and thus leads to an effective decrease of the compression and heating 
factors (l2]i-(|3j as suggested by the SPH results. However the true effect of shock waves was not 
studied in depth and remains up to now still unclear. One could inversely argue that the lack of 
reliability of the SPH models could explain the discrepancy with the AM. 

On the one hand, the simulations had a poor spatial resolution due to the excessively weak 
number of particles used 500 for Bicknell & Gingold [T983l and ^ 7000 for Laguna et al. [T993] 
and Fulbright |1995l l. but also due to an undesirable effect of the smoothing kernels. In the SPH 
approach, the kernels behave as interpolating functions whose size is continually adjusted to keep 
a constant number of particles as interpolation points. In these standard conditions, Fulbright et 
al. (1 19951 1 showed that the usual spherically symmetric kernels were likely to artificially freeze the 
spatial resolution along the direction of compression of the star. To improve this problem, Fulbright 
( 119951 1 resorted to anisotropic spheroidal kernels. Bicknell & Gingold ([1983) also used anisotropic 
ellipsoidal kernels but with too few particles for the kernels geometry to be effective, whereas 
Laguna et al. ( I1993I I used standard spherical kernels. 

On the other hand, the artificial viscosity could produce spurious entropy during the tidal com- 
pression. As is well known, the artificial viscosity is used to numerically deal with shock waves, in 
order to correctly reproduce the shock front without excessive oscillations and smearing, but also in 
order to produce entropy during the propagation of the shock wave. As Fulbright et al. d 19951 1 also 
clearly showed, the usual artificial viscosity was likely to produce entropy during the free fall phase 
even in the physical absence of shock wave development. Fulbright ( 119951 ) implemented a more 
adapted artificial viscosity to reduce this excessive dissipation, but Bicknell & Gingold ( I1983I I and 
Laguna et al. ( 119931 ) used standard forms. Luminet & Carter d 19861 1 already glimpsed this issue 
where including viscosity in the AM, they showed that they could fully reproduce the results of 
Bicknell & Gingold ( 119831 ). The excessively high value of viscosity nevertheless suggested that too 
much dissipation really occured during the SPH simulations. 

Apart from SPH, no other numerical strategy has been set up to follow the three-dimensional 
evolution of a tidally compressed star. Khokhlov et al. ( 119931 ) solved the hydrodynamical equations 
using the Eulerian finite-difference flux corrected transport technique, but they exclusively inter- 
ested in weak penetration factors which did not lead to the compression of the stellar core, and did 
not investigate deeper encounters due to the lack of spatial resolution of the stationary numerical 
grid. Marck et al. d 19961 1 designed a well-adapted formalism combining symmetry considerations, 
an adaptive moving grid, and pseudo-spectral methods, but they exemplified their model to the tidal 
disruption process for a weak penetration factor /? = 1.5, and unfortunately did not performed any 
extensive study of deeper encounters. 

Recently Kobayashi et al. (I2004I I considered a one-dimensional hydrodynamical model based 
on a Lagrangian Godunov-type method (see appendix of their paper) to simulate the direction of 
compression of the star orthogonal to the orbital plane. They clearly showed in the single simulated 
case /3 - 10, through the density and pressure profiles, that a shock wave formed after the pressure 
at the centre of the star reached its maximum. Actually Fulbright ( 119951 ) noticed a quite similar 
behaviour for /? = 5 when the increase of the gas entropy in the stellar envelope occurred after 
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the instant of maximal core compression. Both resuhs then indicate that shock waves are Ukely 
to develop until shortly after the bounce phase as long as y6 < 10. Kobayashi et al. (I2004I I also 
underlined that the propagation of the shock wave within the star could produce a X-ray flare in the 
keV energy, thus leading to interesting observational consequences. 

In this context, the main motivation of the present investigation is to solve the still remaining 
controversy between calculations based on the simplified AM and on the SPH models regarding 
the influence of shock wave development on the compression of the stellar core. We have per- 
formed a detailed hydrodynamical study of the tidal compression process for penetration factors 
3 < j8 < 15 covering all the realistic cases of strong encounters between a main-sequence star and 
a massive BH. We have restricted to follow the evolution of the star from the tidal radius until the 
neighbourhood of the periastron, i.e. during the only free fall and bounce-expansion phases of the 
stellar matter. As justified in the next section, this restriction allows to deal with the tidal compres- 
sion process in a one-dimensional way, considering only the direction of compression orthogonal 
to the orbital plane and neglecting the motion of the stellar matter along others directions. In a 
similar way of Kobayashi et al. ( |2004| l. the hydrodynamical equations have been solved through a 
Godunov-type approach, highly competitive in the numerical treatment of shock waves. Because 
such an approach does not rely on the artificial viscosity technique, the dissipation issue which 
was encountered by the SPH models during the compression of the star cannot arise. Estimates 
of the compression and heating factors at the instant of bounce have been calculated in order to 
investigate the possibility of the thermonuclear explosion of the star. 

In this preliminary study, we have considered minimal assumptions focusing on a Newtonian 
description of the BH gravitational field, and a main-sequence star modeled by a polytropic equa- 
tion of state as was already the case in previous hydrodynamical simulations. The extension to 
general relativistic backgrounds and to a more elaborate equation of state describing the stellar 
structure is postponed to a forthcoming investigation. 

The article is organised as follows. In Sect. 3, we describe the hydrodynamical model and the 
Godunov-type numerical method for dealing with shock waves. In Sect. 4, we present our results 
and show the existence of two regimes depending on whether the shock wave development occurs 
before or after the instant of maximal compression of the stellar core. The critical value of the 
penetration factor separating the two regimes sensibly depends on the assumed equation of state. 
In Sect. 5, we discuss some astrophysical perspectives. 

3. Hydrodynamical model 

3.1. One-dimensional approximation 

If the full tidal disruption process of the star obviously requires three-dimensional calculations, it is 
however possible to reduce the dynamics of the only tidal compression process to one-dimensional 
calculations. 

As previously mentioned, as soon as the star substantially penetrates (say with /? > 3) within 
the tidal radius, the BH tidal field quickly dominates the internal pressure and self-gravity fields. 
The stellar matter gets into a transitional free fall phase in the external gravitational field which 
fully controls its own dynamics. 
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It is well known (see appendix iBli that the tidal field produces stretching along one of the tidal 
principal directions parallel to the orbital plane, and compression at once along the second tidal 
principal direction parallel to the orbital plane and along the tidal principal direction orthogonal to 
the orbital plane (hereafter called the vertical direction). The compressive vertical direction is fixed, 
whereas both compressive and stretching directions within the orbital plane continually change as 
the star moves (see Fig. lB.ll of appendixlB). Moreover, by symmetry relative to the orbital plane, 
the induced vertical motion is fully decoupled from the induced motion in the orbital plane so that 
they can be considered as independent. 

Therefore, moving within the tidal radius until the neighbourhood of the periastron, the star is 
continually compressed along the vertical direction, while in the orbital plane the stretching and 
compression partly cancel out following the rotation through a nearly right angle of both tidal 
principal directions. The AM (Carter & Luminet I1983I I indeed showed that the sections of the 
star in the orbital plane do not significantly change by a factor more than two, which was further 
confirmed by the SPH simulations (Fulbright |1995l ). For these reasons, the dynamics within the 
orbital plane can be quite neglected relative to the dynamics along the vertical direction in a first 
approximation. It has to be noted that during the subsequent bounce phase, when the stellar matter 
vertically expands, the deformations within the orbital plane also remain negligible due to the 
rapidity of the process. 



3.2. Governing equations 

We assume that the motion of the star in the external gravitational field is equivalent to the motion of 
a particle of mass M», since the size of the star is small enough relative to its distance from the BH 
so that the deviation efifect on the motion of its centre of mass is weak. Moreover, the asymptotic 
velocity Voo of the star far from the BH being typically of the order of a few hundred kms 
the specific kinetic energy is much inferior to (the absolute value of) the BH gravitational 
potential at the tidal radius GM,/Rj so that the orbit of the star is assumed parabolic (for a strictly 
parabolic orbit Voo = 0). Note that being only interested in the behaviour of the star inside the tidal 
radius, our results would be applicable to sufficiently elliptic or hyperbolic orbits as well. 

Considering a Cartesian coordinate system (Xi , X2, X3) with origin at the BH and oriented such 
that the motion of the star is contained in the (Xi,X2) plane (Fig.[Tl), the Newtonian parabolic orbital 
law writes 

Xl(t)+l2RlX2(t) = l2Rp(2GM.Rpy'^t, (5) 

X,it)-^-R„ (6) 
X3(t) = 0, (7) 

taking the origin of time f at the periastron of the orbit, and the star-BH distance is 
X{t)^^+Rp. (8) 

The hydrodynamics is described in the reference frame of the star centre of mass. To simulate 
the compression orthogonal to the orbital plane of the star, we consider the vertical direction z 
passing through the centre of mass and joining the poles of the star. 
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Fig. 1. Newtonian parabolic orbit of the star centre of mass S relative to the point-mass BH B. The 
BH is at the origin of the Cartesian coordinate system (Xi , X2, X^) used to follow the star orbit, and 
oriented such that the motion is contained in the {Xi,X2) plane. The distance between the star and 
the BH is X, and the periastron of the orbit is at distance from the BH. The small arrows indicate 
the direction of the star motion as given by (ISll-©. The z axis is attached to the star centre of mass 
and is at each time orthogonal to the orbital plane {Xi,X2). 

The contribution to internal forces is due to the pressure field. The simulation of the tidal 
compression process starts at the tidal radius so that the self-gravity of the star can be neglected 
relative to the tidal contribution. The effects of viscosity, heat conduction, mixing and convection 
within the stellar matter are neglected because their timescales are long enough compared to the 
duration of the compression. 

The contribution to external forces is due to the BH tidal field. The general expression of the 
tidal gravitational acceleration is recalled in appendix lAl From (IA.16b and (IB. Il l, the tidal acceler- 
ation at position z, when the star is at distance X(t) from the BH, is given by 

.fcO = -||z. (9) 

The hydrodynamics is governed by the one-dimensional Euler equations for the conservation 
of mass, momentum, and energy (see e.g. Landau & Lifshitz[l963]l: 

5,p + 5,(pv) = 0, (10) 
d>(pv) + d,(P+pv^)^P8, (11) 
d,(pe) + dA(P+pe)v)^pvg, (12) 

where p is the density, P the pressure, v the velocity, and e = s + the specific total energy 
which includes the specific internal energy s and the specific kinetic energy. 

The stellar matter is considered as a polytropic (ideal) gas of adiabatic index 7: 

P = "RpT (13) 
^pe(y-l), (14) 

where fi is the specific gas constant and T the temperature. The (local) speed of sound is given by 

(15) 

The Euler equations (fT0]i-(fT2ll are expressed in conservation form 




d,U + d:F(U)^S(U,z,t), 



(16) 
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defining 

U=[p,pv,pe]\ (17) 
F(U) = [pv,f + pv^(P + pe)v]^, (18) 
S{U,z,t) = [0,pg,pvg]', (19) 
respectively the vector of conserved variables, of fluxes, and of sources. 

3.3. Numerical method 

The Euler hyperbolic system of conservation laws ( ITSI l is solved by a Godunov-type high-resolution 
shock-capturing finite-volume method (see e.g. LeVeque 120021 Toro ll999l for textbooks). 

Dividing the spatial domain into cells of right (resp. left) boundary z,-t 1/2, centre Zi = (Zi-1/2 + 
Zi+i/2)/2, and width Azi = Zm/i - Zi-iji, 

U(zi,t)=^ r'"'u(z,t)dz (20) 

is defined as the average value over the cell of the conserved variables (fTTI ). and 
lJ"i^V{zi,f) (21) 

as its numerical approximation at time level f" . 

A source splitting approach is used to deal with the source terms, splitting ( fT6] l into the (homo- 
geneous) system of conservation laws 

5,{7-h5,F(1/) = 0, (22) 
and the system of ordinary dififerential equations 

d,U ^S(U,z,t). (23) 

We consider the Strang splitting where the set of numerical solutions |[/"| are advanced from time 
level f" to time level f"^' - f + Af" according to 

Uf' = Xi (Af72) £2 (Af") £1 (Af'72) {t/;'} , (24) 

with Af" the time step constrained by the CFL condition, Xi and X2 the independent solution 
operators for ( l22b and ( l23T l respectively. For clarity, both solution operators are next expressed 
over the time step Af" from time level f" to time level f"^' . 

Numerical scheme for the system of conservation laws 

From the integration of ( l22b over the space-time domain [z,-i/2,Zi+i/2] x [f",f"^'] and using 



U(zi,r')^U{zi,f) + -1 



( rt" 



^ Fmzi-ii2,t))dt 



F{U{zMi2J))dt 



(25) 



The temporal evolution of the numerical solution (1211 1 is then computed by the explicit scheme in 
conservation form 
Af" 

t/«+i ^U"^+— [Pi-,12 - FM12) , (26) 
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where the numerical fluxes 

1 f'"" 

FmH^— y F(U(zM/2,t))dt (27) 

are defined as the approximation of the average value over the time interval [f, f""^'] of the fluxes 
([18]) at the right boundary of the /* cell. 

We have considered the conservation form scheme implemented in the Relativistic Piecewise 
ParaboHc Method (RPPM) of Marti & MuUer ( IT996I . The RPPM actuaUy solves the one- 
dimensional homogeneous Euler equations in special relativity. However, since in our present ap- 
plication of the tidal compression the velocity of the stellar gas is much smaller than the speed of 
light, the relativistic Euler equations reduces to the Newtonian version. Following the Godunov- 
type approach, the numerical fluxes f ,+1/2 are computed in the RPPM by solving at time level t" a 
Riemann problem at the cell boundary Zi+1/2' 

d,U + d,F{U) = 

[/(z,o4^^^'^^ (28) 

with constant states U^^^i^ and f^f^i^,. The RPPM implements an exact solution of the Riemann 
problem (l28l l. For practical reasons, we have instead switched to an approximate solution through 
the use of a Roe-type Riemann solver where the numerical fluxes are given by 



f /.1/2 - \ 



(29) 



- \~^^\h-(u%2-ulu2yr^ 

k=0,± 

with Ak, fic, and respectively the eigenvalues, the right and left eigenvectors (see e.g. Toro |1999] l 
of the Jacobian matrix 

evaluated at the average value U - j (l^J|i/2 + ^i+i/2)- 

Numerical scheme for the system of ordinary differential equations 

From the integration of ( l23T l over the space domain [zi-1/2, Zi+1/2] and using ( |20] |. 

^U(zi,t)^S(zi,t) (31) 
at 



S(zi,t)=— I S(U(zJ),z,t)dz (32) 



where 

1 p+i, 

is defined as the average value over the /* cell of the sources ( fT9] l. 

It is possible to explicitly express 5 as function of U. Let us call Vj the 7* component of any 
vector V. From ( fT9] ) and ( [Tt] ). ( |32] ) writes 



S(zi,t) 



1 



0,— ( "f/i(z,f)^(z,f)* 

T 



U2(z,t)g(z,t)dz 



(33) 
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Since the cell width Az, is small, 

j(z,t)g(z,t)dz 



« — Uj(z,t)dzx— g(z,t)dz (34) 

so that 

S(zi, t) = [O, U,(zi, t) -g{zi, t), U2{zi, t) g{zi, t)^ , (35) 
defining 

g{zi,t) = ^ r'\{z,t)dz (36) 
the average value over the /* cell of the tidal acceleration Q which writes 

= (37) 

with the star-BH distance X{t) evolving with time according to the orbital law Q-®. 

Therefore by ( [3TT l and ([35), the temporal evolution of the numerical solution ( |2TI ) is computed 
from the first-order system of ordinary differential equations 

at 

= [Q,'UMgizi,t),mt)-gizi,t)Y , (38) 

which is solved from time level t" where 1f(t") = U" to time level f"+' where 'Z/(f"+') - 1/"+'. 
We have performed the temporal integration with the explicit fourth-order accurate four-stages 
Runge-Kutta method (see e.g. Lambert ri991| l. 



4. Results 

4.1. Initial and boundary conditions 

In the following, we consider a BH of mass M, - IO^Mq, and a star of mass M» - Mq and 
radius = Rq initially modeled as a polytrope of polytropic index 3/2. The profiles of density and 
pressure along the vertical direction z are obtained from the resolution of the Lane-Emden equation 
(see e.g. Chandrasekhar [T967] l . The profile of temperature deduces from the ideal gas law ( [13] ). 
whereas the stellar gas is uniformly at rest. Since the vertical compression of the star is symmetric 
relative to the orbital plane, the spatial domain is restricted to the positive values < z < 
where the origin corresponds to the centre of the star. We note p, the initial density, P, the initial 
pressure, and the initial temperature at the centre of the star, with p, = 8.43 x lO^kgm"^, 
P. = 8.65 X 10'"* kg m"' s"^, and = 7.53 x lO*" K. The speed of Ught in vacuum is noted c. As 
said before, in order to be able to neglect the self-gravity of the star during the tidal compression 
process, and to study the motion of the stellar matter in a one-dimensional way, the simulations 
start when the star is at the tidal radius ([T|). 

Like in previous SPH calculations, the evolution of the stellar matter obeys the polytropic gas 
equation of state ( [T4l ) of adiabatic index y = 5/3, well-adapted to the structure of small main- 
sequence stars. However we have also interested in cases with y = 4/3 in order to see the influence 
of the adiabatic index on the shock wave development. 
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We consider penetration factors (|4]i of the star within the tidal radius 3 < P < 15. Let us mention 
that the full disruption of a poly tropic star y = 5/3 on a parabolic orbit roughly occurs fory6 > 1. 
In particular, the star is disrupted for fi > 0.7 in the AM (Luminet & Carter [T986] l. whereas it is 
disrupted for p > 0.9 in the improved AM of Ivanov & Novikov ( 120011) . or for /? > 1.5 in the SPH 
model of Fulbright ( |1995l l. Note that the parameter rj = (^Rp/R,^^^^ (M»/M,)'^^ is sometimes used 
to quantify the strength of the star-BH encounter instead of the penetration factor, with the relation 
J] = P^^^^, so that the disruption occurs for rj < \. 

From a numerical point of view, the spatial domain < z < 7?, is divided in sub-intervals 
of unequal length uniformly discretized, and has been covered by «j 8000 cells in the different 
simulations. During the compression of the stellar gas, the more exterior cells of the numerical 
domain become empty, i.e. the density in these cells decreases, the more interior cells fill. A cell 
is considered empty when the density becomes inferior to the critical value lO^^p*, and is taken 
off the numerical domain by imposing vacuum values (finite small values for numerical reasons) to 
the hydrodynamical variables. This condition allows to follow the decrease of the star radius with 
time. Also, ghost cells are added to compute the numerical fluxes ( |29l ) for the cells situated at the 
extremities of the numerical domain. At each time step, reflecting boundary conditions are imposed 
to ghost cells to the left of the numerical domain due to the symmetry relative to the orbital plane, 
whereas vacuum conditions are imposed to ghost cells to the right. 

4.2. Case y = 5/3 

4.2.1. Description of tine tidal compression process 

The hydrodynamical simulations performed for star-BH encounters in the range 3 < yS < 15 have 
underlined that the evolution of the tidal compression process actually depends on the penetration 
factor with respect to a critical value 12 for a polytropic gas y - 5/3. 

Encounters 3 < j8 < 12 

The evolution of the velocity profile for an encounter with y6 = 7 is reproduced on Fig. |2l 
between the instant when the star is initially at the tidal radius and the instant when the compression 
at the centre of the star is maximum. The homologous velocity distribution, i.e. viz, t) ~ z, at the 
different times indicates that the contribution of the internal pressure field is fully negligible relative 
to the contribution of the BH tidal field (see Bicknell & Gingold[l983]l. Therefore, the whole stellar 
matter vertically collapses in free fall in the external gravitational field, actually with both subsonic 
and supersonic velocities. The free fall velocities increase as the star approaches the periastron, and 
leads in the present case to high Mach numbers up to ^ 25 (Fig.|2]bottom). 

The free fafl phase finally breaks up following the sudden increase of the compression at the 
centre of the star (Fig.|3]l, just after the passage through the periastron of the orbit. Indeed the tidal 
field decreases after the periastron, but continues to compress the stellar matter until the central 
pressure suddenly builds up to fully counteract the compressive contribution. The stellar matter 
then bounces and expands with subsonic velocities from the centre, while it continues to collapse 
elsewhere with subsonic and supersonic velocities. The evolution of the hydrodynamical variables 
can be followed on Figs.|4l|5] 
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Fig. 2. Velocity and Mach number profiles along the positive vertical direction z at different times 
t during the free fall phase for/? = 7 and y - 5/3. Labels stand for the following values of t [&]: 
(0) -894.56, (1) -105.02, (2) -58.57, (3) 6.30 x 10"^ (4) 23.19, (5) 44.97. The star is at the tidal 
radius at f = -894.56 and at the periastron at f = 0. The Mach number is defined by |v(z, t)/a{z, t)\, 
with a(z, f) the local speed of sound given by Eq. ( fTSl ). The position Zs of the sonic point where 
k(Zs, 01 - ci(Xs, t) is indicated by the black points. At left (resp. right) to the sonic point, the flow is 
subsonic (resp. supersonic) with |v(z < Zs, f)l < a{z < Zs, t) (resp. |v(z > Zs, 01 > «(z > Zs, t)). The 
homologous velocity profile, i.e. v(z, f) ~ z, shows that the whole stellar matter collapses in free 
fall in the BH gravitational field. 



During that phase, pressure waves are produced and accumulate near the sonic point, where 
they finally steepen into a shock wave (Fig.|6]l. Powered by the expanding motion, the shock wave 
propagates outwards through the collapsing matter, imparts momentum in its own direction (Fig.|4] 
bottom) and heats the medium (Fig. |5]bottom), until it reaches the star radius. The whole stellar 
matter then continues to expand, but we did not further follow the evolution because the one- 
dimensional approximation becomes thereafter less and less valid. 

As summarized in Table [U the formation properties of the shock wave during the bounce- 
expansion phase of the stellar matter actually depend on the penetration factor In particular, as the 
penetration factor increases, the shock wave forms closer to the centre, as indicated by the position 
of the sonic point, and faster after the instant of maximal central pressure. The shock wave then 
propagates within the collapsing matter during a shorter time before reaching the star radius, since 
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Fig. 3. Evolution of the central pressure as function of time t foi'P - 7 and y = 5/3. After the 
passage of the star through the periastron at f = 0, the stellar matter suddenly compresses at the 
centre of the star until the pressure reaches the sufficient magnitude to counteract the compressive 
contribution of the BH tidal field. The central pressure is maximal at f ^ 44.97. 



Table 1. Characteristic quantities when the shock wave forms for penetration factors 3 < yS < 12 
in the case 7 = 5/3. Column 2: time t [s] when the shock wave forms (the star is at the periastron at 
t = 0). Column 3: position of the sonic point Zs [10^^ R,] at time t. Column 4: star radius R [10^^ R^] 
at time t. Column 5: position of the sonic point relative to the star radius. Column 6: time interval 
Afi [s] between the maximum of central pressure and the formation of the shock wave. Column 7: 
time interval At2 [s] between the formation of the shock wave and the end of the collapse of the 
whole stellar matter. 
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the free fall velocities increase with the penetration factor. The evolution of the density and velocity 
profiles in the case /5 - 10 is reproduced on Fig. |7] Let us mention that no shock wave forms for 
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Fig. 4. Density and velocity profiles along the positive vertical direction z at different times t during 
the bounce-expansion phase foryS = 7 and y - 5/3. Labels stand for the following values of t [s]: 
(5) 44.97, (6) 45.51, (7) 46.08, (8) 46.32, (9) 46.56, (10) 46.79, (11) 46.96, (12) 47.03. The star 
is at the periastron at f = 0. The central compression is maximal at f » 44.97. The stellar matter 
expands from the centre of the star while it collapses elsewhere. A shock wave forms in front of 
the expansion and propagates outwards. The collapse stops at f 47.04. 



The propagation velocity of the shock wave and the compression ratios of the shocked mat- 
ter are given in Table |2] The shock wave propagates within the stellar matter with velocities 
X! 1 0^ km s" ' , and produces density jumps tending towards the strong shock limit {y+l)l{y-l) - ^ 
(see e.g. Landau & Lifshitz [l963l l. 

Such an evolution of the tidal compression process until yS < 12 is in full agreement with the 
qualitative results of Fulbright (1 19951 1 and Kobayashi et al. (I2004I I mentioned in Sect. |2] in the 
respective cases [3-5 and [3 - 10. It also confirms for these star-BH encounters the assumption of 
the AM where no shock wave forms before the instant of maximal compression of the stellar core. 

Encounters 12 < j8 < 15 

When the star penetrates within the tidal radius with a penetration factor y6 > 12, the hydrody- 
namical simulations have shown that a shock wave could indeed form during the free fall phase. 
The evolution of the hydrodynamical variables is given on Figs.[8]|9]for an encounter with fi - 12. 
As the whole stellar matter collapses in the BH gravitational field, in contrast with weaker encoun- 
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Fig. 5. Pressure and temperature profiles along the positive vertical direction z at different times / 
during the bounce-expansion phase for j0 = 7 and y = 5/3. Same labels as in Fig.lD 

Table 2. Characteristic quantities of the shock wave for penetration factors P - 1 and 10 in the 
case y - 5/3. Column 2: compression ratio (Xshock (defined as the ratio of density just upstream 
of the shock front to density just downstream of the shock front). Column 3: velocity of the shock 
front Vshock [10^ km s The values of time t [s] correspond to those of Fig. |4]for yS = 7, and of 
Fig.|2]for/3= 10. 
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ters the profiles progressively steepen until they finally break into a shock front. The shock wave 
then propagates up to the centre where it collides with the symmetric shock wave propagating on 
the other side of the orbital plane, which produces an additional compression of the central matter. 
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Fig. 6. Velocity profiles along the positive vertical direction z at different times t during the bounce- 
expansion phase for P - 1 and y = 5/3 (enlargement of Fig. |4]bottom). The stellar matter sub- 
sonicaUy expands from the centre of the star while it both subsonically and supersonically col- 
lapses elsewhere. Both opposite motions produce pressure waves which propagate outwards with 
the speed of sound, and which accumulate near the sonic point (black points) where they steepen 
into a shock wave. 

It is interesting to notice that the shock waves precisely form during the sudden increase of the 
compression at the centre of the star (Fig.fTOt induced by the tidal field after the passage through 
the periastron. 

After collision, both shock waves are reflected outwards which indeed stops the central com- 
pression. The evolution of the hydrodynamical variables during the bounce-expansion phase can be 
followed on Figs. [TTIfT2l The propagation of the shock wave reverts the collapse, while the stellar 
matter expands from the centre under the effect of the pressure. The collapse definitely stops a few 
time after the instant of reflexion when the shock front reaches the star radius. As in the previous 
description, the shock wave speeds up the expansion in the less dense regions (Fig.[TT]bottom) and 
heats the collapsing matter (Fig.[T2lbottom). 

The formation properties of the shock wave as function of the penetration factor are summarized 
in Table [3] The evolution of the density and velocity profiles during the free fall and bounce- 
expansion phases in the case [i - \5 are also reproduced on Figs. [T3]fT4l Actually, the shock wave 
forms further from the centre of the star as the penetration factor increases. It therefore propagates 
on a longer distance before reflexion so that the shock front more increases. We have noticed in 
these different star-BH encounters that the shock wave forms when the central pressure * 10^ - 
10* P. (seeFig.[IO|. 

The propagation velocity before and after reflexion of the shock wave and the compression 
ratios of the shocked matter are given in Table |4] 
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Fig. 7. Density and velocity profiles along the positive vertical direction z at different times t during 
the bounce-expansion phase foryS =10 and y - 5/3. Labels stand for the following values of t [s]: 
(1) 20.80, (2) 20.91, (3) 21.02, (4) 21.14, (5) 21.25, (6) 21.37, (7) 21.41. The star is at the periastron 
at f = 0. The central compression is maximal at f » 20.80. The stellar matter expands from the 
centre of the star while it collapses elsewhere. A shock wave forms in front of the expansion and 
propagates outwards. The collapse stops at f w 21.42. 



It is interesting to notice that the tidal compression of stars presents a qualitative analogy with 
the core collapse of massive stars (i.e. type Ib/c and type II supernovae) in the sense that both pro- 
cesses involve homologous velocity profiles. However they are significant physical differences. 
In the hydrodynamical supernovae scenario, the collapse of the stellar core is due to the self- 
gravitational field and is spherically symmetric. The core splits into a subsonically homologous 
collapsing inner part and a supersonically (non homologous) falling outer part (Yahil & Lattimer 
I1982I I. The shock wave development is due to the stiffening of the equation of state, i.e. the in- 
crease of the adiabatic index from a value somewhat below 4/3 to a value ^ 2.5 - 3, when the 
central density reaches the nuclear matter density, which stops the collapse from the centre of the 
core and produces pressure waves moving outwards and accumulating near the sonic point where 
they steepen into a shock front (see e.g. Muller [r998l l. Moreover, after having bounced behind the 
shock wave, the inner core quickly settles into hydrostatic equilibrium, and the supersonic accre- 
tion continues for many dynamical timescales. On the other hand, in the tidally compressed stars 
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Fig. 8. Density and velocity profiles along the positive vertical direction z at different times t during 
the free fall phase foryS =12 and y - 5/3. Labels stand for the following values of f [s]: (1) 13.40, 
(2) 13.73, (3) 13.85, (4) 13.94, (5) 13.99. The star is at the periastron at f = 0. The position Zs of 
the sonic point where |v(zs, 01 = <2(zs, f), with aiz, f) the local speed of sound given by Eq. (fTsT l. 
is indicated by the black points. At left (resp. right) to the sonic point, the flow is subsonic (resp. 
supersonic) with |v(z < Zs, t)\ < a{z < Zs, (resp. |v(z > z^, t)\ > a{z > Zs, t)). As the whole stellar 
matter collapses in free fall in the BH gravitational field, a shock wave forms and propagates up to 
the centre of the star. 



scenario, the collapse is triggered by the external (tidal) gravitational field of the BH and occurs 
vertically in the direction orthogonal to the orbital plane. The core only collapses homologously, 
whereas the shock wave development is due to the sudden build-up of central pressure. 

4.2.2. Compression and heating factors 

With a view to subsequent consideration on the thermonuclear reactions that can be triggered within 
the tidally compressed star, particular interest attaches to the maximal values of the central density 
and temperature. 

The behaviour of the central density and temperature is reproduced on Fig.[T5]for the star-BH 
encounters previously considered. 
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Fig. 9. Pressure and temperature profiles along the positive vertical direction z at different times / 
during the free fall phase for^S -Yl and j - 5/3. Same labels as in Fig. [8] 

Table 3. Characteristic quantities when the shock wave forms for penetration factors 12 < p < 15 
in the case y = 5/3. Column 2: time f [s] when the shock wave forms (the star is at the periastron 
at f = 0). Column 3: position of the shock front formation Zshock [10^^ R*]- Column 4: star radius 
R [10^^ Rt] at time t. Column 5: position of the shock front formation relative to the star radius. 
Column 6: time interval Afi [s] between the formation of the shock wave and the maximum of 
central pressure. Column 7: time interval At2 [s] between the formation of the shock wave and 
the end of the collapse of the whole stellar matter. Column 8: time interval Atj, [s] between the 
formation of the shock wave and its reflexion at the centre of the star 
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The evolution as function of the penetration factor, of the maxima pf and respectively 
reached by the central density and temperature at the instant of maximal compression, and of the 
duration Af™ during which these values are maintained, is given in Table |5] and is reproduced on 
Fig. [El 



20 



M. Brassart and J. -P. Luminet: Shock waves in tidally compressed stars 



10^ 



^ 10^ 



CO 
CO 

o 



Q. 



CD 



10' 








25 



30 



35 



5 10 15 20 

Time [s] 

Fig. 10. Evolution of the central pressure as function of time f for/3 = 12 and y = 5/3. After 
the passage of the star through the periastron at f = 0, the stellar matter suddenly compresses at 
the centre of the star During that phase, a shock wave forms at f x 13.96 (a) on both sides of the 
orbital plane (see Figs.[8]|9]l. Both symmetric shock waves propagate inwards until they collide at 
the centre of the star at f » 13.98 (b), which produces an additional (instantaneous) compression of 
the central matter (b) —> (c). The shock waves are then reflected outwards which stops the central 
compression at / » 14.10 (c). 



Let us recall that the AM (Luminet & Carter [T986] l predicted for a polytropic gas y = 5/3 that 
the compression and heating factors of the stellar core were respectively given by 
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For comparison, the fit of results of Table|5]to power laws leads to the evolution laws 
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Fig. 11. Density and velocity profiles along the positive vertical direction z at different times f 
during the bounce-expansion phase for/? =12 and y = 5/3. Labels stand for the following values 
off [s]:(6) 14.04, (7) 14.11,(8) 14.18,(9) 14.29,(10) 14.41,(11) 14.46. The star is at the periastron 
at f = 0. The central compression is maximal at f 14.10. The shock wave formed during the free 
fall phase, after reflexion at the centre of the star at f 13.98, propagates outwards which reverts 
the collapse. Behind the shock wave, the stellar matter expands from the centre under the effect of 
the central pressure. The collapse stops at f « 14.47. 

There is actually an excellent agreement between both results ([39l)-(l4T]) and (|43])-(|45]) for the 
whole star-BH encounters 3 < < 15. As recalled in Sect. |2l the AM assumes that the vertical 
collapse occurs without shock wave development until the bounce of the stellar core. We have 
clearly shown that this assumption is indeed verified for star-BH encounters with /? < 12. On 
the other hand, for stronger encounters, the shock wave development as the whole stellar matter 
collapses has no particular effect on the compression and heating factors. This is due to the fact 
that the shock waves form late enough during the sudden central compression, so that when they 
stop the process after collision and reflexion at the centre of the star, the central matter has already 
been compressed enough by the tidal field. 

We have not taken into account the thermonuclear energy generation during the tidal compres- 
sion process in the present study. However, since we have confirmed the core thermodynamical 
quantities predicted by the AM at the instant of bounce, we can be confident that the correspond- 
ing results akeady obtained on the possible thermonuclear reactions remain valuable in the range 
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Fig. 12. Pressure and temperature profiles along the positive vertical direction z at different times 
t during the bounce-expansion phase for /3 - 12 and y - 5/3. Same labels as in Fig.fTTI 



3 < /3 < 15. In this regard, coupling a nuclear network to the AM, Luminet & Pichon (1 19891 1 in- 
terested in the nucleosynthesis in a main-sequence stellar core approaching a BH of 10^ - 10^ M©. 
They calculated that, starting from a typical population I chemical composition (i.e. a mixture of 
hydrogen, helium, carbon, nitrogen, and oxygen), the main nuclear reactions during the bounce 
phase will be accelerated proton captures on seed elements, that will be achieved for^S > 10. An 
interesting point is that nucleosynthesis processes do not occur only during the maximal core com- 
pression when the density and temperature are very high, but will also continue during the cooling 
expansion phase, mainly due to weak decay of unstable isotopes on timescales of few tens of sec- 
onds. The total thermonuclear energy release largely exceeds the gravitational binding energy of 
the star, so that the dynamics of the stellar debris will be modified in the sense that a larger fraction 
of the ejected gas could be unbound to the BH. In a forthcoming paper we shall examine in more 
detail how the hydrodynamical model with shock waves influences the nucleosynthesis processes. 



4.3. Case y = 4/3 

In order to see the influence of the compressibility of the stellar gas on the shock wave development, 
calculations have been also performed for the polytropic gas y - 4/3. This case may give an 
interesting indication on what will happen with a more realistic equation of state. Indeed for the 
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Fig. 13. Density and velocity profiles along the positive vertical direction z at different times t 
during the free fall phase for/? - 15 and y - 5/3. Labels stand for the following values of t [s]: 
(1) 8.13, (2) 8.36, (3) 8.47, (4) 8.59, (5) 8.71, (6) 8.77. The star is at the periastron at f = 0. The 
position Zs of the sonic point where |v(Zs, 01 - ciizs, t), with a(z, t) the local speed of sound given by 
Eq. ( fTsT l. is indicated by the black points. At left (resp. right) to the sonic point, the flow is subsonic 
(resp. supersonic) with \v{z < Zs,OI < ciiz < 2s,0 (resp. |v(z > Zs,f)l > > Zs,t))- As the whole 
stellar matter collapses in free fall in the BH gravitational field, a shock wave forms and propagates 
up to the centre of the star 



strongest star-BH encounters, it is clear that the compressed stellar core will become more and more 
radiative owing to the excessively high values of density and temperature reached. For this reason, 
Luminet & Carter ( 119861 ) already considered an initial mixture of a non-relativistic polytropic gas 
together with a black body photon gas within the AM. They nevertheless found results, in particular 
the compression and heating factors of the stellar core, which were intermediate between those 
obtained with the pure polytropic gases y - 5/3 and y - A 13. 

The hydrodynamical simulations have shown that the tidal compression process occurs in the 
same way as the case y = 5/3, and that there exists a critical value of the penetration factor separat- 
ing the shock wave development before or after the instant of maximal compression. Precisely in 
the case y = 4/3, shock waves form during the bounce-expansion phase of the stellar matter when 
3 < /3 < 6, whereas they form during the preceeding free fall phase when 6 < /3 < 15. 
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Fig. 14. Density and velocity profiles along the positive vertical direction z at different times f 
during the bounce-expansion phase for/? =15 and y = 5/3. Labels stand for the following values 
of t [s]: (7) 8.82, (8) 8.89, (9) 8.94, (10) 8.98, (11) 9.02, (12) 9.03. The star is at the periastron at 
f = 0. The central compression is maximal at f ^ 8.81. The shock wave formed during the free fall 
phase, after reflexion at the centre of the star at t ^ 8.80, propagates outwards which reverts the 
collapse. Behind the shock wave, the stellar matter expands from the centre under the effect of the 
central pressure. The collapse stops att^ 9.04. 



The evolution of the compression and heating factors as function of the penetration factor are 
given in Table |6] and are reproduced on Fig.[T71 The fit of results to power laws correctly leads to 
the evolution laws 
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Table 4. Characteristic quantities of the shock wave for penetration factors (3 - 12 and 15 in the 
case y - 5/3. Column 2: compression ratio cTshock (defined as the ratio of density just upstream 
of the shock front to density just downstream of the shock front). Column 3: velocity of the shock 
front Vshock [10"' kms"']. The values of time f [s] correspond to those of Figs. [8] and [TTI for fi - \2, 
and of Figs.[T3]and[T3for/3 = 15. 
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Table 5. Evolution of thermodynamical quantities at the centre of the star as function of the pen- 
etration factor 3 < y6 < 15 in the case 7 = 5/3. Column 2: time t [s] when the central density and 
temperature reach their maximal value (the star is at the periastron at f = 0). Column 3: maximal 
value pf [p»] of the central density. Column 4: maximal value {Tt\ of the central temperature. 
Column 5: duration Af" [lO^'^r.], with t. = (1/Gp*)'^^ x 1.33 x 10^ s, of the maximal central 
density and temperature. It is defined as the full width at half maximum of the peak of central 
temperature. Column 6: star radius R [10^^ R,] at time t. Column 7: density pR [lO^'^p,] at the star 
radius. Column 8: temperature Tr [10^^ r»] at the star radius. 
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Fig. 15. Evolutions of the central density and temperature as function of time f for different pen- 
etration factors p in the case y = 5/3. After the passage of the star through the periastron at 
f = 0, the stellar matter suddenly compresses at the centre of the star (see Fig. [3] for P - 1 and 
Fig. [To] for y6 = 12). The increase of central density and temperature reaches its maximum at 
t ~ 44.97,20.80, 14. 10, 8.81 for respectively /? = 7, 10, 12, 15. 




The AM (Luminet & Carter [T986] l however predicted for the whole star-BH encounters that 
— 1 ^ 0.007 /3^ (50) 

P* ) AM 

^ * 0.19/?2. (51) 

I ^ • J AM 

In contrast with the poly tropic gas 7 = 5/3 where the hydrodynamical and AM results are quite 
consistent, the shock wave development during the free fall phase (i.e. for /3 > 6) with y - 4/3 has 
for effect to substantially decrease the central density and to slightly increase the central temper- 
ature for the strongest encounters. In such cases, it appears that shock waves indeed form earlier 
during the sudden central compression (Figs. n"8lfT9] l, so that the central matter is less compressed 
by the tidal field when the shock waves stop the process after reflexion outwards from the centre 
of the star. The amount of internal energy of compression lost is however balanced by a slightly 
superior amount of energy carried by the shock waves propagating inwards until reflexion. 
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Fig. 16. Evolutions of the maximal central density, maximal central temperature, and duration dur- 
ing which the maximal values are maintained as function of the penetration factor y6 in the case 
y - 5/3. Triangle symbols: values computed by the hydrodynamical model (see Table|5]l. Dashed 
line: power law predicted by the AM (see Eqs. ([39]l-(l4TI)). In bottom figure, the characteristic inter- 
nal timescale of the star t, = (1/Gp,)'^^ ~ 1.33 x 10^ s. 



5. Summary and perspectives 

When a star deeply penetrates within the tidal radius of a massive BH without falling in the event 
horizon, the tidal effects produce a short-lived strong compression of the stellar core, at the end 
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Table 6. Evolution of thermodynamical quantities at the centre of the star as function of the pen- 
etration factor 3 < jS < 15 in the case y = 4/3. Column 2: time t [s] when the central density and 
temperature reach their maximal value (the star is at the periastron at ? = 0). Column 3: maximal 
value [pt] of the central density. Column 4: maximal value Tf [T*] of the central temperature. 
Colunm 5: duration Af™ [IO^'^t.], with t» = HIGp*)^^^ ~ 1.33 x 10^ s, of the maximal central 
density and temperature. It is defined as the full width at half maximum of the peak of central 
temperature. 
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of which it becomes highly flattened according to a pancake-shape configuration. The first in- 
vestigations performed in the 1980's within the framework of the semi-analytical AM predicted 
that the compression could actually trigger some thermonuclear reactions within the steUar core 
which were indeed likely to lead to an explosion of the star. Afterwards, three-dimensional hydro- 
dynamical simulations based on SPH carried out in the 1980-1990's came to the conclusion that 
the pancake stars phenomenon could be less spectacular than originally believed, finding lower 
estimates of the pancake quantities as to the compression and heating factors of the stellar core, 
and therefore questioned the explosive tidal disruption. The authors of these studies claimed that 
the development of shock waves during the tidal compression were mainly responsible for the dis- 
crepancy. However, the SPH results remained uncertain owing to potential issues of the numerical 
method to accurately simulate the compression phase, but also owing to a poor spatial resolution. 

In this article, we have performed a detailed hydrodynamical study of the vertical tidal compres- 
sion undergone by the star in the direction orthogonal to its orbital plane, through a one-dimensional 
hydrodynamical model based on the high-resolution shock-capturing Godunov-type approach to 
correctly deal with shock waves. The main results are the following. 

Firstiy, we have estabhshed a more reaUstic scenario of the tidal compression process, and 
put in evidence two regimes depending on the strength of the star-BH encounter defined by the 
penetration factor 3 < /? < 15. In the case of main-sequence stars evolving through a polytropic 
gas equation of state of adiabatic index 5/3 and for encounters such as /3 < 12, the stellar matter 
collapses in free fall towards the orbital plane until the central pressure suddenly increases to revert 
the process. The stellar matter then bounces at the centre of the star before expanding. Shock waves 
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Fig. 17. Evolutions of the maximal central density and temperature as function of the penetration 
factor y6 in the case 7 = 4/3. Triangle symbols: values computed by the hydrodynamical model (see 
Table|6]l. Dashed line; power law predicted by the AM (see Eqs. (l50ll-(l5Tli). 



develop during the expansion phase on both sides of the orbital plane and propagate outwards. The 
description of the free fall phase until the instant of maximal compression is indeed in agreement 
with the initial assumption of the AM. On the other hand, for stronger encounters, shock waves 
develop during the free fall phase on both sides of the orbital plane, while the central pressure 
suddenly increases, and propagate inwards until they collide at the centre of the star. The shock 
waves are then reflected and propagate outwards, which reverts the free fall phase and causes the 
stellar matter to expand from the centre under the effect of the pressure. 

Secondly, estimates of the compression and heating factors in both regimes are the same as 
those predicted by the semi-analytical model, thus showing that shock waves have no particular 
influence on the compression of the stellar core contrary to what concluded in the SPH studies. 
Therefore we fully confirm the possible triggering of thermonuclear reactions within the stellar 
core, like they were already investigated on the basis of the semi-analytical model. 

Thirdly, in both regimes of compression the shock waves carry outwards a brief (< 0.1 s) but 
very high (above lOOkeV) peak of temperature from the centre to the surface of the star (see 
Table Q and Fig. |20] i. and could provide an eflicient mechanism to transfer the pancake energy 
into hard electromagnetic radiation. This last result is very promising since it could likely give 
rise to a new type of X- or y-ray bursts, as akeady pointed out by Carter (|1992| l on the basis of 
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Fig. 18. Density and velocity profiles along the positive vertical direction z at different times t 
during the free fall phase for/? - 12 and y - 4/3. Labels stand for the following values of t [s]: 
(1) 9.76, (2) 10.46, (3) 11.63, (4) 12.78, (5) 13.93. The star is at the periastron at f = 0. The position 
Zs of the sonic point where |v(zs, 01 - ci{Zs, t), with a{z, f) the local speed of sound given by Eq. (fTSl) . 
is indicated by the black points. At left (resp. right) to the sonic point, the flow is subsonic (resp. 
supersonic) with |v(z < Zs, 01 < <^(z < Zs, (resp. |v(z > Zs, 01 > «(z > ^s, 0)- As the whole stellar 
matter collapses in free fall in the BH gravitational field, a shock wave forms and propagates up to 
the centre of the star. Note the size of the shock front relative to the case y = 5/3 (see Fig.|8]l. 



qualitative arguments and recently considered by Kobayashi et al. (120041 ) through hydrodynamical 
calculations. 

The present numerical simulations are a first step towards the resolution of the full problem of 
strong star-BH encounters. They have been performed in the case of a main-sequence star orbiting 
a massive BH along Newtonian orbits. A more realistic situation should deal with the relativistic 
gravitational field of a static (Schwarzschild) or a rotating (Kerr) BH, for which it has already been 
shown that a deeply plunging star within the tidal radius undergoes several successive compres- 
sion phases. This will require new calculations taking into account such relativistic effects, and one 
can expect that new hydrodynamical phenomena wiU arise due to the interaction of ingoing and 
outgoing shock waves. Including the self-gravity is also necessary in order to start the simulations 
well outside the tidal radius, and to better describe the post-pancake phase. Moreover, since our 
results confirm the possibility of thermonuclear explosion of the star, it will be necessary to couple 
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Fig. 19. Evolution of the central pressure as function of time f for/3 = 12 and y - 4/3. After 
the passage of the star through the periastron at f = 0, the stellar matter suddenly compresses at 
the centre of the star During that phase, a shock wave forms at f 9.85 (a) on both sides of the 
orbital plane (see Fig.fTsTl. Both symmetric shock waves propagate inwards until they collide at the 
centre of the star at f 14.01 (b), which produces an additional (instantaneous) compression of 
the central matter (b) (c). The shock waves are then reflected outwards which stops the central 
compression at f » 14.04 (c). Note that relative to the case j - 5/3 (see Fig.fTOli. the shock waves 
produce a superior compression at the instant of collision, because they form further to the centre 
(see Figs.lSlandfTSll so that the shock fronts more increase before reaching the centre. 

a nuclear network to the hydrodynamics. The cases of more sophisticated equations of state for 
the stellar core (e.g. including radiative corrections) should be investigated, as well as the inclu- 
sion of non conservative terms (due to thermonuclear energy generation) in the energy equation. 
Eventually, in order to make more reliable predictions which could be compared to the available 
observations of high energy flares in galactic cores, the radiative transfer of energy through the 
disrupted star should be taken into account, as well as the aftermath of the stellar debris. 

The subject of the tidal disruption of stars is called to develop in the near future. Indeed, the X- 
UV flares recently observed in the core of non-active galaxies can be interpreted as a consequence 
of such a phenomenon, thus providing new observational constraints for the models. A number of 
future (low- and high- energy) X-ray all-sky surveys are planned, the first ones to be launched in a 
few years such as XEUS and Constellation X missions, so predictions about the rate of these events 
and signals in the X- and y-ray regimes are of strong interest. By providing a quick localization 
of such events, followed by the detection of the corresponding afterglows in the optical, infrared 
and radio bands, such missions could bring to the understanding of stellar disruptions as much as 
the Beppo-Sax and Swift telescopes have brought to the comprehension of y-ray bursts. Some of 
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Table 7. Evolution of characteristic quantities when the shock front reaches the star radius as 
function of the penetration factor 3 < < 15 in the case j - 5/3. Column 2: time f [s] when the 
shock front reaches the star radius (the star is at the periastron at f = 0). Column 3: star radius 
R [10^^/?,] at time f. It corresponds to the minimal value of the star radius reached during the 
compression. Column 4: temperature Tr at the star radius. It corresponds to the maximal 
value of the temperature at the shock front. Column 5: velocity vr [10^ kms '] at the star radius. It 
corresponds to the velocity of expansion when the shock front reaches the star radius. Note that no 
shock wave forms for^S = 3. 
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the y-ray bursts already observed by Swift and other telescopes could be tentatively interpreted as 
pancake stars. 

A massive BH imbedded in a dense cluster of stars strongly disturbs their dynamics within the 
gravitational influence radius (Frank & Rees ll976] l 
GM. 



. ,^ig / M. \/200kms- , 

3 X 10'^ cm I . , (52) 



where Vk, is the stellar velocity dispersion in the galactic core. Currently, it is widely accepted the 
existence of an observational-based law connecting M, and Vm (Ferrarese & Ford 12005]) such that 

If the star is diff'used along an orbit within the tidal radius 
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it is ultimately disrupted by the tidal gravitational field. However, since the BH gravitational radius 
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Fig. 20. Evolution of the temperature as function of time f foryS = 7 and y = 5/3. Ttie star is at 
ttie periastron at f = 0. The solid line corresponds to the temperature at the centre of the star The 
dashed line corresponds to the increase of temperature produced by the shock wave as it propagates 
outwards within the stellar matter (see Fig. |5]bottom). The duration at half maximum of the peak 
of temperature at the shock front is w 0.05 s. 



grows faster than the tidal radius with the BH mass, the tidal disruption can only occur outside 
the horizon for BH mass below the Hills limit (Hills [1975] ). which is a 10^ Mq for a star of solar 
mass and radius. The different characteristic distances are reproduced on Fig. [21] Our calculations 
require BHs in the range 10^ - 10^ Mq, i.e. typically the massive BHs expected in ordinary galactic 
nuclei, such as our own Galactic Center. 

Since the cross section for penetration of the star inside distance Rp is 



the rate of stars penetrating by a factor /3 within the tidal radius is (Luminet & Barbuv lTgQOI l 
dN(/3) 



dt 
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where A^c is the number density of stars at core radius of the stellar cluster As mentioned in Sect.[Tl 
the coefficient 4 x 10 could be considerably increased by the presence of a self-gravitating ac- 
cretion disc around the BH. The important thing is that the frequency of strong star-BH encounters 
(say with > 3) is by no means negligible, just one order of magnitude less than the frequency 
of tidal disruptions. Thus the rate of pancake stars should be about 10"^ per galaxy per year If the 
associated flares are generated in hard X- or y-ray band, the full observable universe is transparent 
and several events of this kind per year would be expected. 
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log(M.) 

Fig. 21. Characteristic distances (in centimeters) for solar-type star / massive BH (in solar mass 
units) interactions. is the gravitational influence radius (|53] l, Rg is the gravitational radius ( |55] ). 
Rt is the tidal radius (|54] |. Also represented is R^, the radius below which the star penetrates within 
Rj by a factor yS > 10 and is likely to detonate thermonuclear reactions. The domain of astrophysical 
relevance for such a process is the grey area. 
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Appendix A: Derivation of tlie tidal gravitational acceleration 

We recall the general expression of the tidal gravitational acceleration from which ^ is derived. 
As in the main text of the article, the BH gravitational field is considered in Newtonian dynamics. 

Let us introduce (see Fig. lA.U a first Cartesian frame (B,exi,ex2,ex3) with origin B at the 
BH position (reference frame S of the BH), and a second Cartesian frame (S ,eri,er2,er3) with 
origin S at the star centre of mass (reference frame S of the star centre of mass, which relative to 
the inertial reference frame S is parallely propagated along the star motion). At time f in reference 
frame S, the star centre of mass is at position X(t) - Xi{t)exi, and a point-mass of the star at 
position x(t) = Xi(t)exi, where hereafter the summation convention over repeated indices is used 
and indices vary from 1 to 3. In reference frame S, the same stellar point-mass is at position 
rit) = nit) en with r,(f) = x,(f) - Xi(t). 

Let us place in reference frame S, and let us compute the (gravitational) acceleration at position 
x{t) in the star due to the instantaneous action of the BH gravitational field. At distance x = (x,x,)'^^ 
from the BH, the stellar point-mass is exposed to the Newtonian gravitational potential 

<^M^-^. (A.1) 
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Fig.A.l. Cartesian reference frame of the BH {B,exi,ex2,ex3) and Cartesian reference frame of 
the star centre of mass (S ,eri,er2,er3). The star centre of mass is located by the position vector 
X(t) = Xi(t)eXi. A point-mass of the star is located by the position vector x(t) = Xi(t)eXi relative 
to the BH, and by the position vector r(t) = rj{t)eri relative to the star centre of mass with r(t) — 
x{t)-X{t). 

Since in the present case, the star radius is very small compared to the distance X = 
between the star centre of mass and the BH, r ^ X with r = {ririyi^. Therefore, the gravitational 
potential (lA.ll l can be expanded on the neighbourhood of the star centre of mass X: 



^{x) = <b{X + r) 

= ^{X) + {xi-xd[d,,<i>{x)]^^x 

+ - Xdixj - Xj) [d,d,^^{x)\^^ 

+ hxi - Xi){xj - Xj){xk - Xk) [d,_d,d,,^{x)\^^ 



+ 0{\xt). 

The expansion iA.2\ is rewritten in the form 

<^{x) = OCX) + {xi - X>) [d,,^(x)]^^^ + OtW, 

defining the tidal gravitational potential 

(Dt(a:) = - ^(xi - Xi)(xj - Xj)CijiX) 

- - Xdixj - Xj)(xk - Xk)Dijk(X) 
+ 0(\xf), 

along with the second-order tidal gravitational tensor 

C,v(Z)^-[5,,5.,<DWL, 

and the third-order deviation gravitational tensor 

A,,(Z) = -[5,,5,,5,,<5(^)L^^. 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
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The deviation tensor (IA.6I 1 measures, at the lowest order, the deviation of the motion of the star 
centre of mass from the motion that it would have if all the mass of the star was concentrated at 
this point. One can indeed show that the equation of motion of the star centre of mass is rigorously 
given by 

^ = -MAMX) + ^Dijk J rjn p d'r + 0(\rf) (A.7) 
with 

dxMX) = ^Xi, (A.8) 

and where the integration is performed over the volume of the star of density p. However since 
Rt <K X, the star appears as a point-mass relative to the BH so that the deviation effect is weak 
enough for the deviation tensor to be completely neglected. The motion of the star centre of mass 
can therefore be considered as strictly equivalent to the motion of a particle of mass M, in free fall 
in the external gravitational potential 'l'(Z), and ( IA.7l i then limits to the first term given by ( IA.8I ). 
In the same way, the tidal potential (IA.4b limits to the tidal tensor term (IA.5b . 

The first and second-order partial derivatives of the gravitational potential ( lA.lb are equal to 

d,,<:>{x) = ^Xi, (A.9) 
GM. ( 3xiXj\ 

5,5,<DW = -^[5,-—^j, (A.10) 
where dij denotes the usual Kronecker delta: 

\o i* J. 

Using dA.ll ), ( |A.9b , and dA.lOb applied at the star centre of mass X, one can express with ( |A.4b 
and ( IA.5I ) the gravitational potential ( IA.3I ) at position x. From that result, relative to reference frame 
S the acceleration at position x deduces directly: 

^ = -5, OCX) (A. 11) 

= - [d.,^{x)]^^^ + {xj-Xj)Cij{X) (A. 12) 

with 



[dx,^ix)]^=x^^^i^ (A. 13) 



GM . 

and 



GM. I 3XiXj 



Cu(X)-^[-Su+^j- (A. 14) 

The first term of ( IA.12b given by (IA.13I ) is identical to the acceleration ( lA.Sb of the star centre 
of mass in the BH gravitational field. The tidal gravitational effects, which are by nature dififerential 
effects, are given by the second term of ( IA.12b referred to as the tidal gravitational acceleration. 
Therefore in the reference frame S of the star centre of mass, the acceleration at position r = x - X 
is equal to the tidal acceleration: 
d^rj _ (fxi^ _ d^Xj 

dt^ ~ dt^ dt^ ^ ■ 

= rjCij(X) (A. 16) 

using ( |A.12b with ( |A.13b . and jATTj l with ( IATSI i. 
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Following the assumptions mentioned in appendix lAl we recall the general properties of the tidal 
gravitational field. As made explicit by the expression ( IA.I6I 1 of the tidal acceleration, the behaviour 
of the tidal field is imposed, actually at the lowest order, by the tidal tensor ( IA.14I ). 

The expression of the tidal tensor can be simplified by choosing the orientation of the coordinate 
axes {exi,ex2,ex3) such that the motion of the star centre of mass remains in the (exi,ex2) plane 
during the encounter with the BH. This choice, which is of course possible since the Newtonian 
orbit of the star is planar, implies that XT,it) = at each time t. The expression ( |A.14| i of the tidal 
tensor then gives 



C(Z) = 



GM. 



-1 + 



X2 





X2 



-1 + 



3X2 
X^ 











-1 



(B.l) 



Moreover, in order to easily show the induced effects by the tidal field on the star, let us place 
in the principal directions frame of the tidal tensor. The latter being of real values, and by definition 
symmetric, can always be diagonalized. The eigenvalues of ( IB. lb are 
2GM. 



X3 ' 

GM. 



^^3 - ^2, 

and the normalized (right) eigenvectors associated with the eigenvalues ( IB.2b -( lBr4l i are 



IX2I (Xi 
eui = — —exi + ex2 

X \A2 



eu2 



X 

eu3 = ex3. 



Xi 

exi + ex2 

X\ 



In the principal directions frame (S,eui,eu2, euj), the tidal tensor is diagonal: 



Cij 




From ( IA.16b the tidal acceleration at position u(t) - Ui{t) eui in the star is equal to 

(f'Ui 

which with ( IB. 81 ) simplifies to 

(f'Ui 



= '^j «y ^ij 



(B.2) 

(B.3) 
(B.4) 

(B.5) 

(B.6) 
(B.7) 

(B.8) 

(B.9) 
(B.IO) 



From (IB. 10b . it is directly noticeable that one of the effects of the tidal field is to stretch the 
stellar matter along the direction eu^, since by ( |B.2b the eigenvalue is positive. On the contrary, 
the eigenvalues A2 and given by ( IB. 31 ) and ( IB. 4b being negative, a second effect of the tidal 
field is to compress the stellar matter along both directions eu2 and eu^,. Actually, as can be seen 
from ( IB. 5b and ( IB.6I ). both principal directions eui and eu2 are contained in the plane of the star 
orbit {exi,ex2)- Because these directions depend on the star centre of mass coordinates X,(f), they 
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continually change with time f (see Fig. IB. It . It is different for the third principal direction eu^, 
which from ( IB.7b is orthogonal to the orbital plane, but also constant so that this direction remains 
always fixed during the motion of the star. 




-100 -80 -60 -40 -20 20 40 60 80 100 



Fig. B.l. Representation of the eigenvectors ( IB.SI l and (IB.6I 1 of the tidal gravitational tensor (IB. Il l 
at different positions along the star orbit. The eigenvectors eui and eui indicate the principal di- 
rections of the tidal tensor and are contained in the plane of the star orbit, eui (dashed arrow) is a 
stretching principal direction whereas eu2 (solid arrow) is a compressive principal direction, and 
both continually change as the star moves. The orbit of the star centre of mass (solid line), whose 
coordinates relative to the BH (cross) are {X\ , X2), is assumed parabolic and is computed with Q- 
(|6l). In this example, M. = IO^Mq, M, - M^, and = Rq. The dot-dashed circle represents the 
BH tidal radius i54\ and the solid circle the BH gravitational radius ( |55] l. The star enters in the tidal 
radius from the right bottom corner with a penetration factor /3 - 5. 
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